function ww = initialize_w(N,K,MG,hGG,pow_max)
% initialize variable
% mW
ww = zeros(MG,N,K);
for nn=1:N
    for kk=1:K
        ww(:,nn,kk) = hGG{nn,nn,kk} / norm(hGG{nn,nn,kk},"fro") * pow_max(nn)/K;
    end
end
end
